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Abstract: I argue that the system of interacting bosons at zero temperature and in ran- 
dom external potential possesses a simple critical point which describes the proliferation 
of disorder-induced defects in the superfluid ground state, and which is located at weak 
disorder close to and above one dimension. This makes it possible to address the critical 
behavior at the superfluid-Bose glass transition in dirty boson systems by expanding around 
the lower critical dimension of d = 1. Within the formulated renormalization procedure near 
d = 1 the dynamical critical exponent is obtained exactly and the correlation length critical 
exponent is calculated as a Laurent series in the parameter ^/e, with e = d — 1: z = d, 



V = l/y/3e + 0{^/e) for the short-range, and z = 1, u = \j2/3e + 0{^/e), for the long-range 
Coulomb interaction between bosons. The identified critical point should be stable against 
the residual perturbations in the effective action for the superfluid, at least in dimensions 
1 < d <2, for both short-range and Coulomb interactions. For the superfiuid-Mott insulator 
transition in the system in a periodic potential and at a commensurate density of bosons I 
find u = {l/2^/e) + 1/4 + 0{^/e), which yields a result reasonable close to the known XY 
critical exponent in d = 2. The critical behavior of the superfluid density, phonon velocity 
and the compressibility in the system with short-range interactions is discussed. 

PACS: 74.40+k, 5.30Jp, 5.70Jk 



1 Introduction 



The transition between the superfluid (SF) and the locahzed insulating (Bose-glass, BG) 
ground states in the disorder bosonic systems at zero temperature is a paradigmatic example 
of a quantum phase transition 0. It is believed to be relevant for a number of experi- 
mentally observed superconductor-insulator transitions at low temperatures, including those 
in Josephson junction arrays 0, amorphous superconducting films 0], underdoped high- Tc 
cuprates p , '^He and ^He in disordered media , [0 , and possibly in two-dimensional elec- 
tron gas in Si-MOSFETs The vortex transitions in type-II superconductors in magnetic 
field in presence of correlated disorder have also been related to this problem 0, [|l^. De- 



spite of its ubiquity in nature and good qualitative understanding of the insulating phase 
0, the critical behavior at the SF-BG transition is only poorly known. This may appear 
surprising, particularly when contrasted with the greatly successful theory of classical critical 



phenomena |]rT| , even in the presence of quenched disorder [|I^ . The source of the difficulty 
can be traced to the difference between the classical and the quantum disordered systems 
that becomes most apparent in the path integral formulation of the quantum problem: d- 
dimensional quantum system at T = maps onto a d+l-dimensional classical theory, but the 
disorder potential appears random only in (i-dimensions, while it is completely correlated in 
the remaining dimension which represents the imaginary time. This is because static disor- 
der is represented by a time-independent random potential. The unpleasant consequence of 
this fact is that there is no upper critical dimension for the SF-BG transition in the standard 
sense. One way to try to deal with this problem is to demote a number of imaginary-time 
dimensions from unity to a small number e,-, and then to perform a double expansion in both 
er and e = A — d — Er |13]. Unfortunately, the perturbative results for the critical exponents 



obtained this way appear to have a very poor convergence in e,-, rendering this approach of 
limited use for calculations of the critical exponents. That this may not be a coincidence 
could be suspected on the basis of certain exact inequalities which should be satisfied by the 
exponents at the disordered critical points ||I4|, and which dictate that their values should 
be very different than at the Gaussian fixed point. For example, m d = 2 and e,- = 1, at 
the SF-BG critical point one expects > 1, and z ^ 2, in contrast to z/ = 1/2 and z = 1 at 
the Gaussian fixed point for the relativistic bosons. In fact, within the double e-expansion 
the disordered and the Gaussian fixed points generally do not coalesce at the upper critical 
dimension dup = 4, but instead the disordered critical point stays always non-trivial (except 
when e = e,- = 0), and becomes unstable only in dimensions higher than some dc > dup [|15|. 
All this is of course completely different from the standard 4 — e expansion for the classical 
systems where the Wilson-Fisher critical point coincides with the Gaussian fixed point 
at the upper critical dimension, becomes unstable right above it, and the critical exponents 
in the physical dimensions are numerically not too different from their Gaussian values. One 
might speculate that the failure of the double e-expansion may be due to non-analyticity of 
the critical exponents as functions of e,- |T^. 



Recognizing the limitations of the double e-expansion many authors turned to different 
ways of addressing the SF-BG transition: large-N |T^, real-space methods strong- 
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coupling expansion [0, numerical calculations and the renormalization group in fixed 
dimension |^ have all been applied to this problem. In absence of a small parameter one is 
left feeling uneasy with the results obtained by these methods, and it would still be desirable 
to have a controlled analytical approach to the SF-BG critical point which would be free of 
the conceptual and the calculational problems of the double e-expansion. Besides providing 
a practical tool for calculation of the universal quantities at the SF-BG phase transition, 
such a scheme would be expected to be able to correctly discern different universality classes 
of superconductor-insulator transitions, similarly to the celebrated e-expansion in the field 
of classical critical phenomena. 

In this paper I present an attempt to develop such a systematic theory of the SF-BG 
critical point as a perturbative expansion in deviation of spatial dimensionality d of the 
quantum system from unity. That the parameter e = d — 1 effectively controls the value of 
disorder at the superconductor-insulator critical point has been suggested previously by the 
author in where this was used to study the superconductor- Anderson insulator transi- 
tion in the model of disordered attracting spinless fermions near one dimension. The basic 
idea is rather general and may be understood qualitatively as follows. Consider a quantum 
system of interacting bosons (which may be Cooper pair composites, for instance) in d = 1 
and at T = 0: in absence of any external potential the Bose single-particle operator will 
exhibit a power-law quasi-long-range correlations, with a non-universal power K depending 
on the microscopic interactions. K may therefore be understood as an exactly marginal 
coupling constant. What will be the effect of a random external potential? It must suppress 
superfluidity by trying to localize the elementary excitations. In the language of renormal- 
ization group this means that it should produce a positive term in the recursion relation 
for in attempt to drive the Bose single-particle correlation function shorter ranged. The 
renormalization group flow of the disorder coupling constant D will in turn be affected by 
so that if K is larger than some K* , small D is relevant and vice versa |^, |^. The 
crucial point is that the marginality oi K m d = 1 without disorder forces the SF-BG fixed 
point to lay at = (and some finite K = K*, see Figure l.b ) in d = 1. To use this 
fact to develop an expansion in small parameter e = d — 1 for the problem with disorder 
one must also observe the following: power-law superfiuid correlations (without disorder) 



are characteristic of c? = 1 [|25[; in d > 1, true long range order will develop at T = 0. In 
d = 1 + e one may still think of K, but no longer as a marginal coupling; it should slowly 
(controlled by e) scale to zero, so that the single-particle correlation function tends to a finite 
constant at large distances. Thus a small deviation of dimensionality of the system from 
unity should produce a small negative term in the recursion relation for K, which then may 
be balanced by the positive term induced by disorder. This opens a possibility of having a 
random critical point in the physical region of the coupling space which is infinitesimal in 
disorder slightly above one dimension (see Figure 1). 

The central assumption made here is that the SF-BG critical point in d = 2 smoothly 
evolved from the trivial (in the sense that D* = 0) SF-BG critical point in d = 1. The 
additional motivation for this idea comes from the study of the disordered Bose-Hubbard 



model [21 1 . Using the duality transformations to map the dirty-boson model on lattice onto 
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a field theory for defect degrees of freedom, it was noticed that the SF-BG critical point in 
the dual formulation arises through similar renormalization scenarios in ci = 1 and d = 2. 



The form of the dual theories in one and two dimensions is quite different [21|, reflecting 



more complicated nature of defects in higher dimensions, nevertheless, in both cases it is 
the couphng of the external potential to the field that describes topological defects of the 
phase of the superfiuid order parameter which ultimately determines the critical behavior. 
The physical mechanism for destruction of superfluidity at T = is postulated to be the 
same in d = 1 and d = 2: the unbinding of phase slip-antislip (in d=l) or vortex- antivortex 
(in d = 2) pairs induced by disorder destroys the phase coherence. The critical point which 
describes this transition is strongly coupled in ci = 2, but becomes weakly coupled in d = 1, 
which suggests an expansion around the latter solvable case. 

In the rest of the paper this idea is used to calculate the critical exponents at the SF- 
BG transition in the systems of disordered bosons with short-range and long-range Coulomb 
repulsion near one dimension. My starting point is the effective low-energy action for the 
phonon excitations above the SF ground state [^] in presence of an external random poten- 
tial. Near d = 1 the non-trivial part of the recursion relations is determined from the theory 
as precisely in d = 1, while the effect of dimensionality enters through the non- vanishing 
canonical dimensions of the coupling constants. To handle the disorder term I take advan- 
tage of the density representations of the theory |^ in d = 1, where it may be understood 
as dual to the phase representation ||2l|, and in which the external potential couples to the 



operator that creates 2n phase slips in the superfiuid ground state. 1 first study the case 
of the short-range interaction between bosons in detail, later to extend the calculation onto 
the Coulomb universality class. The recursion relations are derived to O(e^) and the cor- 
relation length critical exponents u and the anomalous dimension t] are calculated as series 
in e to two lowest orders, whereas the dynamical critical exponent is obtained exactly. In 
particular, within the presented formalism the conjecture of Fisher et al. [0], that the 
dynamical exponent z = d and z = 1 for the short-range and the Coulomb repulsion between 
bosons, respectively, holds exactly as long as the critical point is smoothly connected to the 
one in d = 1. For the short-range interactions, both the superfiuid density and the phonon 
velocity renormalize to zero at the SF-BG transition, as power-laws with powers that vanish 
as ~ ^/e close to d = 1. Compressibility on the other hand approaches a finite constant at 
the transition, in accordance with the expected gapless nature of the Bose glass phase. 

The paper is organized as follows. In the next section I define the problem and determine 
the canonical dimensions of the coupling constants in the effective action in the case of 
short-range interactions. In section 3 the recursion relations for the coupling constants to 
two lowest orders in e = d—1 are presented, and the critical exponents for the case of short- 
range repulsion are obtained. In section 4 I examine the relevance of the residual terms in 
the superfiuid effective action at the SF-BG critical point. Section 5 presents the results 
for the critical exponents in the case of Coulomb repulsion between the bosons. I test the 
expansion method against the known results for the XY universality class relevant for the 
superfiuid-Mott insulator transition in a periodic potential in section 6. Summary and the 
discussion of the main points is given in the last section, and the calculational details are 



3 



relegated to the appendices. 



2 Hamiltonian and the canonical dimensions 

I begin by defining the system of interacting bosons with the standard second-quantized 
Hamiltonian: 

H=^J rf'x|VvI/(f)|2 + 1 1 d'xd''yv{S-m{^)\'\'^{y)\', (1) 

where v{r) is a short-range repulsive interaction, with v{q = 0) finite, = and 
satisfies standard bosonic commutation relations. The Hamiltonian (1) has a superfluid 
ground state in all > 1. The form of the effective theory for the low-energy excitations 
may be derived by writing the Bose operator in the density-phase representation as \l/(a;) = 
p^^'^{x) eyipi(f)(x) ||26[, where p{x) = Po + H(x), with H(x) being a local deviation from the 
average density of bosons po- Expanding the Hamiltonian to the quadratic order in H and 
(j), the imaginary time quantum mechanical action S = J dfixdriji^'^dr'^ + /7[\E'^, at 
finite temperature T becomes: 

S4,,n = j c/'^x^^cir(^(V0(f,r))2 + 27^^;;v^2(f,r) + ^H(x,r)9.0(f,r) + ...), (2) 

where 13 = h/ksT, and the terms of higher order in H(x) and </)(x) and their derivatives 
have been omitted for the moment. The boundary conditions satisfied by the fluctuating 
fields are (f){x,P) = 4>{x,0) + 2'7in{x), with n{x) integer, and H(x, /3) = H(a;, 0). At T = 0, 
the action in Eq. (2) may be understood as an effective low-energy theory for the superfluid 
by defining the coupling constants as vj = hTrps/m and vn = l/(47r/iK), where ps is the 
superfluid density and k is the compressibility of the system. The effective theory S^p^u then 
represents an exact description of the low-momentum {\k\ < A, where A << 27rpo is an 
ultraviolet cutoff) phonons, which are the Goldstone modes of the broken U{1) symmetry in 
the superfluid ground state. It is just the Landau quantum hydrodynamics for the superfluid 
'^He [^. The presence of the last imaginary term in Eq. (2) implies that in the subspace of 
low-energy states the density and the phase should be considered as canonically conjugate 
variables: 

[Uix),<l>iy)] = -t6'ix-y). (3) 

The terms of higher order in H and omitted in Eq. (2) describe the residual interactions 
between the phonons and the deviation of phonon dispersion from linearity. They will be 
considered in Section 4, where I will argue that they are irrelevant at the SF-BG critical point, 
at least for 1 < < 2. Note that these perturbations lead to infrared-finite corrections to the 
low-energy theory for the superfluid, which I will assume are already included in definitions 
of the coupling constants Vj and Vn- 

Let me determine the canonical dimensions of the coupling constants vj and vn from 
the effective theory (2) for the non-interacting excitations at T = 0. Anticipating the effects 
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of disorder, consider the change of these couphngs under rescahng of lengths, but allowing 
for the dimension of the imaginary time to be [r] = — z (that is r ~ L^, where L is a length), 
with z as an undetermined dynamical exponent. The canonical dimensions are 

[vj] = z + d-2, [vn] =z-d, (4) 

as follows from recalling that (f){x) is a phase, so [0] = 0, and [H] = d. The reader recognizes 
the first equation as essentially the Josephson scaling relation for the superfluid density, 
Ps ~ ^"i-d-z gjj^ while the second one is the analogous relation for the compressibility, 
K ~ ^^"'^ 0]. ^ is the diverging correlation length at the transition. Note that the canonical 
dimensions of both couplings vanish in d=l, since for the non- interacting excitations z = 1. 
The last statement is equivalent to the algebraic decay of the single-particle correlation 
function in ci = 1 and at T = at large distances, or, more precisely 

(^+(f)*(o)) ~ poi^r^"^/"'^'^' (5) 

when \x\ —>■ oo. In d > 1 the single-particle correlation function in Eq. (5) tends to a finite 
constant at large separations. 

Consider adding a term with an external potential to the Hamiltonian (1): 

Hr = J d'^xV{x)\^{x)\'^, (6) 

where V{x) is a random function of the coordinate. Evidently, the external potential couples 
directly only to the particle density and not to the phase. One may still perform the Gaussian 
integration over 11 in the action (2), to obtain a complex effective action for the phase 0, 
with disorder entering via an awkward imaginary term i J V{x)dr4>{x,T). Since disorder is 
static, after integration over imaginary time one finds that this gives a purely boundary 
contribution, in which the external potential couples to the windings {n{x)) of the phase in 
the imaginary time. Fortunately, in dimensions d = 1 and d = 2 there exists a real (dual) 
representation of the effective action in terms of the particle density [^, [Q, in which 



disorder appears in a less forbidding form. Here I follow Haldane , and in d = 1 introduce 
a new field 6'{x, r) in the action (2) as n(x, r) = 7i~^dxO'{x, r). Performing the average over 
disorder using replicas after integration over the superfluid phase the effective replicated 
low-energy theory in d = 1 becomes 

N N 

Se = KY, dx dT{c\dx9if+{dr9iY)-DplJ2 dxdrdr' cos2{9i{x,T)-ej{x,r')), (7) 

i=l ^ i,j=i 

where 

9,(x,T) = ei{x,T) + -^ r V{z)dz, (8) 

and I introduced the standard combinations of the coupling constants K = 1/{2ttvj), 
= AvnVj for later convenience, c represents the velocity of the phonon excitations, and 
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K is inversely proportional to the superfluid density. The limit — at the end of cal- 
culations is assumed. Haldane's heuristic derivation of the effective action Sq for bosons in 



one dimension |2i], [27] is presented in Appendix A. It is worth mentioning that the action 
in Eq. (7) also arises as the long-distance theory of the disordered Bose-Hubbard model in 
one dimension ||21|, |20|. Disorder produces an effective interaction between replicas, which 
in d = 1 is reminiscent of the interaction term in the sine- Gordon theory. There are two 
important differences however. First, the interaction is non-local in imaginary time, which is 
a consequence of the quantum nature of the problem, and breaks the relativistic invariance 
of the pure system. Second, the interaction term is invariant under a shift of all fields by an 
arbitrary function of the coordinate x. Both features will have important consequences for 
the critical behavior at the SF-BG phase transition. 



3 Recursion relations and the critical exponents 

In the rest of the paper I will be interested in the Bose system at zero temperature, and 
therefore set /3 = oo in Eq. (7). Without disorder, the effective action in Eq. (2) describes 
the non-interacting excitations. Under the change of the cutoff A — A/s, or equivalently 
under the rescaling of coordinates as x ^ sx and r — > s^r, coupling constants K and c in 
general scale as determined by their relations to vj and vn and their canonical dimensions 
in Eqs. (4), with z = 1. Disorder introduces effective interactions between the excitations 
in any dimension, and in d = 1 the most relevant one has a simple form given by Eq. (7). 
Also, precisely in c? = 1, both K and c have their canonical dimensions vanishing. Thus 
the strategy will be to calculate the recursion relations for K and c in = 1 perturbatively 
in disorder from the action (7), and then to account for the effects of dimensionality when 
d > 1 hj adding terms with the canonical dimensions and with a general z. The reader will 
note that this is precisely the logic of the minimal substraction scheme when used together 
with dimensional regularization in the theory of classical critical phenomena. Finally, I will 
determine the value of the dynamical exponent by requiring the existence of the fixed point 
of the recursion relations in d > 1. 

I now proceed with the implementation of this program, and in d = 1 integrate out the 
modes of the field 6{x,t) which have the momentum A > \k\ > A/s, for ln(s) infinitesimal, 
and with any imaginary frequency — cxd < a; < oo. The details of the perturbative calculation 
to two lowest orders in disorder are presented in Appendix B. Including the d-dependent 
canonical dimensions from the previous section and redefining the couplings as ttK K, 
ttDpq/A^ ^ D, I find the recursion relations: 

{2-z-d)K + ^^+0{D% (9) 



dln{s) ' ' Kc^ 

[z-l-4^)c + 0{D\ (10) 



d\n{s) ' ^2^4 
'^^ {d + 2z- ^)D + 0{D^). (11) 



d\n{s) Kc' 
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The last three equations represent the central result of this paper, and some remarks may 
be in order. To the second order in disorder D the coefficient in front of {drOY in Eq- 
(7) becomes renormalized by the integration over the fast modes, whereas the one in front 
of {dxOY , which is the inverse of compressibility, does not. The recursion relation for the 
velocity c is thus completely determined by the one for K, to two leading orders. To complete 
the scaling transformation I need to specify the dynamical exponent: I choose z to keep the 
velocity c constant under renormalization. This is essential to find a fixed point in the theory 
and also conveniently decouples the recursion relation for c from the other two. Numerical 
value of c is then completely arbitrary. Setting c = 1 and inserting z = 1 + AD / K"^ + 0{D'^) 
into Eqs. (9) and (11), I find the fixed point of the above renormalization group equations 
at 

^* = ^-| + 0(62), (12) 

D* = ^ + 0{e'). (13) 

At the fixed point the dynamical exponent thus equals 

z = l + e + 0(e3). (14) 

Note that there is no O(e^) correction in the last result. The reason is the absence of 
renormalization of the coefficient in front of {d^OY-teTm in the action (7) to the second order 
in disorder. In Appendix C I show that the exact symmetry of the disorder term in the one- 
dimensional theory in Eq. (7) under 9i{x, r) — > Oi{x, r) + f{x) for arbitrary f{x), guarantees 
that this coefficient remains unrenormalized to a// orders in D in d = 1. Within the proposed 
calculational scheme where the non-trivial part of the recursion relations follows from the 
theory as in c? = 1, this represents the sufficient condition for that all higher corrections to 
z in Eq. (14) vanish, i. e. 

z = d (15) 

exactly. The conjecture of Fisher et al. thus follows as an exact consequence of the 
formulated renormalization procedure near one dimension. 

Linearizing the renormalization group flow around the critical point in Eqs. (12)- (13) 
(Figure 1) one finds the correlation length exponent which determines the flow along the 
relevant direction to be 

u= ^ + Oi^/~e). (16) 



/3e 

The other eigenvalue which determines the flow along the SF-BG separatrix gives the leading 
correction to scaling close to d = 1. u turns out to be a Laurent series in y/e and therefore 
of weaker dependence on dimensionality than one would expect. For e = 1 the result (16) 
yields u ~ 0.58, smaller than in the Monte Carlo calculations of Wallin et al. pO[, where 



u = 0.9 ± 0.1. The generalized Harris criterion [jl4| which requires that u > 2/d at & 
disordered fixed point is also not satisfied at this order. It is possible that the higher order 
terms in the expansion will increase the value of z/, and remedy both discrepancies. The 
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divergence of as e — > is a sign of vicinity of the lower critical dimension for the SF-BG 
transition. 

Right at the critical point in c? > 1 the average single-particle correlation function shows 
an algebraic decay at large distances (see Appendix D): 



G(f) -pol^r^*'*- (17) 
The anomalous dimension rj is defined by K*c* = d + z — 2 + 7] [0], so that: 

V = l-j + 0{e'). (18) 

For e = 1 the above result yields a negative r], as one would expect based on the assumption 
that the single-particle density of states iV(e) oc e ^ at the critical point diverges at zero 
energy 0. The result however is much more negative than the one obtained by Wallin et al 



in their Monte Carlo calculation i] = —0.1 ± 0.15. Higher order terms may bring these 
two estimates closer, although expansions of the anomalous dimension around the lower 
critical dimensions are usually notoriously badly convergent. 

The reader may have noticed that the coefficient in front of the second term in the 
scaling equation for K depends on the precise definition of the coupling D and is therefore 
arbitrary. It is reassuring to realize that none of the exponents in fact depends on the value 
of this coefficient. 

Upon approaching the SF-BG transition from the superfluid side, the velocity of phonons 
and the superfluid density in ci = 1 + e for e small both vanish as 

c~5v^, (19) 

and 



Ps ~ (20) 

where 5 is a parameter that tunes through the transition. In contrast, ind = 1 both quantities 
approach flnite values at the superfluid side of the transition, and then vanish discontinuously 
in the insulating phase. Compressibility k ~ Ps/c^, on the other hand, remains flnite at the 
transition in all dimensions, as implied by the found exact relation z = d. Although flnite, at 
the transition its value remains non-universal and dependent on the microscopic couplings 
in the system. 



4 Residual terms in the effective action 

For the flxed point in Eqs. (12)-(13) to govern the critical behavior at the SF-BG transition 
ind > 1 one must assure that the terms of higher order in Il{x) and (j){x) and their derivatives 
which were dropped in the effective action S^^u are irrelevant perturbations. These terms 
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originate from the kinetic energy part of the Hamiltonian (1), which in the density-phase 
representation takes the form: 



r2 

Hk — 



The effective action for the superfluid will thus contain additional parts coming from the 
second and the third term in the above expression: 

-^^.n = / d''xdT[v,n{x){V<l>{x)f + Vy{VU{x)Y + 0{n{VU)% (22) 

with Vx and Vy as two new coupling constants, which together with S'^^n in Eq. (2) determine 
the full effective low-energy action for the elementary excitations, without disorder. 

To examine the relevance of the two new terms in S'^ ji at the SF-BG critical point in 
d > 1 we first need the canonical dimensions of Vx and Vy from Eq. (22): [vx] = z — 2 and 
[vy] = z — 2 — d. As long as z = d at the SF-BG critical point, a small Vy which describes 
the deviation of the spectrum from linearity will stay irrelevant close to d = 1, where its 
dimension may change form the canonical one only by e. It thus seems reasonable that 
this will remain so in physical dimensions d — 2 and d — 3. The coupling constant Vx on 
the other hand, under the assumption that z = d at the SF-BG fixed point is irrelevant by 
power-counting close to d = 1, but becomes marginal in d = 2 and relevant above d = 2. The 
situation is analogous to the effect of the term in the Ginzburg-Landau- Wilson theory 
at the Gaussian fixed point, and one must examine the corrections to the dimension of Vx 
introduced by disorder. Dropping the I'j^-term in Eq. (22), after integration over the phase 
in d — 1 one finds a term 

s', = AnvxK' J dxdT{drey{{dxe) + oiidxef)}, (23) 

in addition to So in Eq. (7). To the two lowest orders in disorder I found previously that 
only the coupling K in front of {drOy-term becomes renormalized. Under the change of 
cutoff A ^ A/s then 

VxK' ^ Vx{s)K'{s) = VxK' + 0{D\ vl). (24) 
Inserting K[s) I then find the recursion relation for Vx'- 

dvx ^ IQD 



dln{s) ^ K^d 



+ 0{D'))vx + 0{vi). (25) 



When d — 2 and z — 2, the fixed point D* /K*'^ ~ 1, so to truly answer the question 
of relevance of small Vx at the SF-BG critical point we need all the higher order terms in 
D in the bracket in the equation above. However, the first term in the expansion of the 
scaling dimension of Vx in D suggests that Vx most likely stays irrelevant even for d = 2. 
The negative sign of this term comes from the effective increase of K under the infiuence 
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of disorder, which on physical grounds may be expected to persist to higher orders. It is 
plausible then that is irrelevant at the SF-BG critical point, at least for 1 < d < 2. 

Finally, consider the higher harmonics in the effective action in Eq. (7) that arise from 
the components of the random potential with the wave-vectors k ~ 2m7rpo with m > 2 (see 
Appendix A): 

oo ^ r 

SeM = - ^ J2 / dxdrdr' cos2m{9i{x,T) - 6j{x,T')). (26) 
Completely analogously as for D, one finds the recursion relations: 

"^^^ {d + 2z-^)D^ + 0{Dl), (27) 



d\n{s) Kc 

so that all Dm with m > 2 are irrelevant at the SF-BG critical point, similarly to the 
irrelevance of vortices with higher vorticity at the Kosterlitz-Thouless transition of the 2D 
XY model. 



5 Long-range Coulomb interaction 

Now that the formalism is in place it is interesting to consider the case of long-range Coulomb 
interactions between bosons in the Hamiltonian (1), Vc{r) = e^/r. The quadratic part of the 
long-distance action at T = is now 

SsTi= I d'^xdT{ — {V(t){x,T)f + / d'^x' —^—Uix,T)Uix\T) + m{x,T)dMx,T)) (28) 
J 271 J \x — x\ 

for d > 1. The canonical dimension of coupling vj thus remains the same as in the case of 
short-range interaction, while the canonical dimension of the charge e can be read from the 
Eq. (28): 

[e^] =z-l. (29) 
I define the long-range interaction in c/ > 1 as 

Precisely in ci = 1 this definition coincides with the short-range and not with the Coulomb 
interaction, which has the Fourier transform ~ ln(l//c) m d = 1. Since I am not interested 
m d = 1 per se, but only in c? > 1, for the present purposes this distinction is not important 
||33|| . The only difference in the scaling equations from the short-range case is then in the 
replacement of the recursion relation (10) with 

> + ^^--f7T-l)c + 0{D^), (31) 



ciln(s) ' 2 K'^c^ 
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where K is defined exactly as before, but = Avje^. Tlie Eq. (31) may be understood 
as tlie recursion relation for charge, in place of the equation for the phonon velocity in the 
short-range case. Assuming that at the SF-BG fixed point the value of the charge is finite 
p8| fixes the value of by requiring c = 1 for example. Combining Eqs. (31), (9) and (11), 



similarly to the short-range case I find the dynamical exponent now to be 

z = l. (32) 

This is again an exact result within the present formalism, for the same reason as in the 
short-range case. The correlation length exponent is 

+0{^e). (33) 



V3e 

For e = 1 this yields v ~ 0.82, within the bounds set by the Monte Carlo calculations of 
Wallin et al. v = 0.90 ±0.15. The inequality u > 2/dis not satisfied by the lowest order 



estimate, but presumably higher-order terms would correct this. The anomalous dimension 
of the single-particle propagator is now 

r/ = ^-^ + 0(6^). (34) 

The result for e = 1 is again much smaller than the Monte Carlo estimate |^ where 
rj = 0.8±0.4. It may be interesting to note that the difference between anomalous dimensions 
for Coulomb and short-range interactions to the lowest order in e is roughly the same as found 
numerically. 

The analysis of the relevance of the residual interactions in the superfiuid action from the 
preceeding section can easily be repeated in the present case with the same final conclusion. 
In fact, since one expects that z = 1 for the Coulomb interaction, both and Vy appear 
irrelevant in d = 2 and d = 3 already by simple power counting. 



6 Superfluid-Mott insulator transition in periodic po- 
tential 

It is well known that in a purely periodic potential, Bose system at integer number of bosons 
per period undergoes a transition from a superfiuid into an incompressible Mott insulator 
(MI) with increasing repulsion between bosons, which is in the universality class of the XY 
model ||3^. Since the critical exponents of the XY model are known virtually exactly ||35[ 



one 



may test the presented method of calculation on this standard case. In ci = 1 the effective 
action in Haldane representation becomes 

Sg = 2nT^y f dxdT{{dJ{x,T)f + {dr9{x,T)f}-v [ dxdr cos{2n{po-l)x + 29 {x,t)), (35) 
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so that for one boson per period (po = 1) it reduces to the sine-Gordon theory. Relativistic 
invariance of the theory imphes immediately that z = 1. In d = 1 + e the renormahza- 
tion group equations for the SF-MI problem are simple deformations of the celebrated 
Kosterlitz-Thouless recursion relations 



"^'^ -e(l -u) + 2v^ + Buv^ + 0{v\ mV) (36) 



(iln(s) 

dv 



=2uv + Av^ + 0(uv^,ev), (37) 
aln(s) 

where u = (1 — l/SnTxy), T^y is the temperature of the equivalent classical XY model in 
d = 2 and v is the fugacity of vortices. The coefficients in Eqs. (36) and (37) are known to 
be i? = 4 and A = —5 [Q, but I left general A and B in the recursion relations to emphasize 
that these two quantities, in contrast to the coefficients of the lowest order terms, separately 
are not universal. Linearization of the flow around the simple critical point of the recursion 
relations (36)-(37) then yields the correlation length exponent 

^.y = ^(1 - ^^^^^v^+ 0(e)). (38) 

It is interesting that the exponent u to this order depends only on the universal linear 
combination of the coefficients A and B ||3^. This is a non-trivial check of the consistency 
of the expansion around the lower critical dimension for this problem, and suggests that the 
same universality may be expected in the SF-BG case as well. Taking e = 1 in the last 
equation one finds u^y = 0.75, reasonably close to the 3D XY value of 0.67 For e = 2, 
Uxy = 0.60, somewhat farther from the expected Gaussian value 1/2. 



7 Summary and discussion 

In this paper I have presented a new analytical approach to the critical behavior at the 
superfluid-Bose glass transition at T = based on the expansion around the lower critical 
dimension d = 1. In the vicinity of d = 1 the SF-BG critical point becomes infinitesimal in 
disorder, which enables one to obtain the recursion relations for the coupling constants in 
the theory as expansions in small parameter e = d — 1. 1 found that the critical exponent 
is z = d {z = 1) for the short-range (Coulomb) interaction between bosons, in accord with 
the conjectures based on finiteness of the compressibility (charge [^]) at the transition. I 
obtained the correlation length exponent and anomalous dimension to two leading orders for 
the short-range and the Coulomb universality classes, and demonstrated that the analogous 
calculation for the superfluid-Mott insulator transition leads to the results in reasonable 
agreement with the known u^y exponent in d = 2 + 1. Higher-order anharmonic terms in 
the effective action for the superfluid which describe interactions between the excitations 
and the deviation of the dispersion relation from linearity are argued to be irrelevant at the 
SF-BG critical point, at least for 1 < d < 2. 
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The critical point found in this work should be understood as describing the proliferation 
of topological defects in the superfiuid ground state, which are induced by the disorder 
potential at T = 0. In this sense the SF-BG (and even more SF-MI) transition is similar to 
Kosterlitz-Thouless finite temperature transition in the XY model; both are defect-mediated 
transitions, except that the role of entropy term in the quantum (T = 0) problem is played 
by the external (random or periodic) potential. This is precisely why the dual formulation 
of the problem is useful: it rewrites the theory in terms of " correct" degrees of freedom. 

A similar expansion around one dimension for the superconductor-insulator transition 
has been attempted earlier by Kolomeisky who argued that the physical critical point 



exists only in < 1, in contradiction with my conclusions. The mistake in his work lies in 
using the effective action in c? = 1 in Haldane form (7) to infer the canonical dimensions of 
the coupling constants in d ^ 1. The transformation from the particle density to the variable 



6' is a particular realization of the duality transformations ||2^, |31], and in this form 
it is possible only in ci = 1 and meaningless everywhere else. For this reason one needs to 
determine the canonical dimensions of K and c from the effective action in the form (2), 
which is general and exists in any dimension. Note that this does not contradict my later 
use of the effective action in Eq. (7) to calculate the recursion relations, since the non-trivial 
parts of these are determined by the theory precisely in d = 1. 

The reader should note that the considered SF-BG transition in dirty bosonic system 
should be the right universality class for the s-wave superconductor - localized insulator 
transition in disordered electronic systems, at least in vicinity of one dimension. This is 



know to be true precisely in ci = 1 [^, and should remain true in its neighborhood as well. 



since the attractive interaction in the singlet channel will still open the gap in the spin sector 



Experimentally, the simplest universal quantity to measure is the product of dynamical 
and correlation length exponents, which follows for example from collapsing the finite tem- 
perature resistivity data on either insulating or superconducting side onto a universal curve 
Although the experimental situation is presently somewhat unclear, most measurements 
are consistent with z ^ 2 and > 1 for the short-range interactions between bosons 0], [Q. 
It appears that before attempting a more realistic comparison with experiments the present 
calculation would need to be pushed to higher order. This may be done by remaining faithful 
to the logic of minimal substraction scheme which would require the higher order terms in 
the recursion relations to be calculated only in d = 1. The results for the SF-MI transition 
in the periodic potential, where I relied on the already calculated higher order recursion 
relations, are encouraging in this sense. 

The presented theory also facilitates a systematic treatment of number of interesting 



issues in the field, like the interplay between periodic potential and disorder |^ , |^ , the fi- 
nite temperature crossovers close to the SF-BG criticality [^], and the corrections to scaling, 
which are left for future considerations. 

It is a pleasure to acknowledge useful discussions with I. Affleck, M. Sharifzadeh Amin, 
V. Barzykin, P. Stamp, Z. Tesanovic, M. Oshikawa and A. Zagoskin. The author is grateful 
to NSERC of Canada and the Izaak Walton Killam foundation for financial support and to 
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Aspen center for physics where this work was initiated. 



8 Appendix A: derivation of the effective action in d = 

1 

Haldane proposed a way to improve upon the hydrodynamic theory in Eq. (2) in d = 1 
by including the structure of the quantum fluid at shorter length scales. Here I outline 
the reasoning which leads to the effective action (7), which forms the basis for most of the 
calculations in this paper. Introduce a field 9"{x) so that dxO"{x) = npi^x). As the density is 
discreet and therefore a sum of delta-functions, 9"{x) jumps by n at the positions of particles. 
So we may write 

oo 

p{x) = n-^dj"ix) 5[9"ix)-nn], (39) 

n=— oo 

or using the Poisson summation formula and introducing the local deviation from the average 
density 11 (a;): 

oo 

p(x) = (po + n(x)) E (40) 

m=— oo 

The m = term in the last formula may then be understood as the coarse-grained density, 
and the rest of the sum yields shorter and shorter length-scale corrections. The imaginary 
time action in Eq. (2) together with the random potential term now becomes 

Ssu= [ dx ['^ dT{^{dx(j){x,T)Y + 27rvNll'^{x,T) + iU{x,T)dr^{x,T) + (41) 

J Jo 271 

V{x)n{x, r) + poV^(a;)(e''^"("'") + c.c.)) + hh 

where "hh" denotes the higher harmonics that arise from the infinite sum in Eq. (40). 
Introducing a new field 9'{x,t) as TT~^dx9'{x,T) = II{x,t), one may write the term that 
couples the particle density and the superfluid phase as 



7T 



J dx dT{dx^{x,T)){dr9'{x,T)). (42) 



The phase can now be straightforwardly integrated out which leads to 

Se= Idx f dT(^-^{dJ\x,T)f + -^{dr9\x,r)f + (43) 

J JO TT 27rf J 

V{x)-dx9'{x, t) + po^(x)(e2*(^'("'")+"«'") + c.c.)). 

TT 

Fourier components of V{x) with the wave- vectors k << po couple only to Il{x) in Eq. 
(43), and can be completely absorbed by a shift of 11, i. e. by defining a yet another field 
9 as in Eq. (8). These components lead only to forward scattering, and in c? = 1 are 
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unrelated to localization. The components of V{x) with k ~ 27rpo however couple directly 
to expi2^'(a;), which is the operator that creates 2n phase-slips in the ground state, and can 
not be eliminated this way. Assuming the Gaussian probability distribution 

V*(k)V(k) 

P[V{k)] ~ e ^ (44) 

for the Fourier components with k ~ 27rpo, after introducing replicas and averaging over 
V{k) one obtains the Eq. (7) as in the text. 

For a different, and maybe somewhat more rigorous derivation of the theory (7) the 



reader should consult the second paper in the reference |21 



9 Appendix B: momentum shell RG 

Define the slow and the fast components of the field 6: 

eix,T) = eiix,T) + e2ix,T), (45) 

where 

e,{x,T) = f r ^e{k,co)e^'^^^-\ (46) 

J~A/s 2tT J~oo ZTT 

and analogously the 62, but with the integral over momenta going over A/s < \k\ < A. To 
the lowest order in D, after integrating out the fast modes the remaining action for the slow 
modes becomes 

Se, = KY, J dxdTic\dJi,{x,T)Y + {^r6^,{x,T)Y) - (47) 

i 

DplJ2 I dxdrdr' cos2{9iAx,r) - 9,^Ax,r')){cos2{e2Ax,r) - e2jix,r')))2 + 0{D''), 

where the average is performed over the quadratic part of the action for the fast modes. 
First, consider the terms off-diagonal in replica indices. One finds 

{cos2{e2,i{x,T) - e2jix,r')))2 = s'^^ (48) 

for i j. If i = j 

{cos2{e2M r) - 62^ r')))2 = exp -(^^ H , (1 - e'-(— ')) + 0{Hsf)). 

TV J\ J~oo C li. -\- UJ 

(49) 

The first term in the bracket in the 0(ln(s))-term matches the contribution from the off- 
diagonal terms in Eq. (48) and together they renormalize D: 

D D{s) = /}s(3-(-i^'=)-^) + o(D'). (50) 
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After expanding the cosine-term in Sg^ to the quadratic order we see that the time-dependent 
term in Eq. (49) apart from generating some irrelevant terms also renormalizes the coefficient 
in front of {drOiy. Thus, for an infinitesimal ln(s) 

K K{s) + ln(.) + 0(D^ Hsf). (51) 

To the first order in D the non-local term in x does not get corrected by the elimination of 
the slow modes: 

Xc' ^ K{s)c{sf = Kc^ + 0{D^). (52) 

After redefinitions of the couphng as described right above the Eq. (9) we see that dif- 
ferentiating the above expressions yields the lowest-order terms in the Eqs. (9)-(ll) when 

z = l. 

Next, consider the contribution quadratic in D to the Se^^: 

E / dxdTdT'dydvdv'{cOS2{ei,i{x,T) - e^j{x,T'))cOS2{9,,rn{y,v) - 9,,n{y,v'))^^^^ 



i,j,m,n 



[(cos2(6'2,i(a;,r) - 92j{x,t')) cos2{92,m{y,v) - 92,n{y,v')))2 - 

{cOs2{924x,t) - 92j{x,r')))2{cOs2{92,Uy,v) - 92,n{y,v')))2\ + 

sin2(6'i,j(x,r) - 9ij{x,t')) sm2{9i.^,{y,v) - 9i^n{y,v')) 
{sm2{92,i{x,T) - 92j{x,T'))sm2{92,m{y,v) - 92,n{y,v')))2}. 



Using the identities: 



(cos a cos (5)2 = e 2 cosh(Q;/3)2, (54) 

(sinQ;sin/3)2 = e~ 2 sinh(Q;/3)2, (55) 

and realizing that all averages on the right hand side in the last two equations are ~ ln(s), 
we see that the first term in Eq. (53) is of order (ln(s))^. To the first order in ln(s) only 
the second term in Eq. (53) can potentially contribute to the action for the remaining slow 
modes, and it equals to 

- E / dxdrdr'dydvdv' sin 2{9i^i{x, r) - 9ij{x, r')) (56) 

Sin 2(^i,,(y, v) - 9^,r.{y, ^'))(^2,i(x, r)92,i{y, v))2 + 0{{Hs)f). 

Since the corelator appearing in the last equation drops exponentially in |r — and \x — y\, 
the dominant contribution comes from the leading term in the expansion of the fields under 
the second sinus around the point (x, r) and (x, r'). For n — j one finds the term 

8i:>Voln(s) 



J dxdrdr' sin^ 2 (^1,^(2;, r) - 9ij{x, r')), (57) 



which obviously, besides adding a constant to the free energy, just generates a higher har- 
monic. The recursion relation in Eq. (11) is therefore correct to 0{D^'). There may be 0{D^) 
rcnormalization of K deriving from the term (56), but the reader may convince himself that 
this does not matter for the exponents u and rj to the two lowest orders in e. 
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10 Appendix C: proof that z = d to all orders 



I derive the general form of the recursion relation for the coupling constant c, from which it 
will follow that for the short-range interactions z — d within the considered renormalization 
procedure. First step is to show that the coefficient in front of the {d^OY- term in the theory 

(7) in d = 1, K(? — 2i'jv/vr, can not become renormalized during the integration over the 
fast modes. Let me assume that the Bose system is of a finite size L and impose the periodic 
boundary conditions in space direction: 

^,(0,t) = ^,(L,t) (58) 

for all < T < 1/T. To be precise 1 will assume a finite temperature T, although this does 
not affect the argument in any way. The renormalized vn{L) at the scale L can be defined 
as _ 

where S is the imposed twist in the boundary conditions in the space direction: 

9imr)^9i{0,r)+S. (60) 

Defining a new field 9[{x^t) = 6i{x,T) — 6x/L which satisfies the periodic boundary condi- 
tions, and noticing that the disorder term in Eq. (7) is invariant under this transformation, 
it readily follows that 

vn{L) = vn, (61) 
independently of the length scale L. In other words, under the scale transformation x ^ sx 

d _ 2 



d\n{s] 



Kc' = (62) 



exactly, in ci = 1. If one defines a function F{D, K, c) which is not explicitly dependent on 
dimension as 

dK 

__ = ^2-z-d)K + F{D, K, c) (63) 
in d > 1, then Eq. (62) implies that 

dc , , F(D,K,c), 

Demanding that the last two renormalization group equations vanish at some non-trivial 
couplings D*, K* and c*, implies that at the fixed point 

z^d. (65) 
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11 Appendix D: Bose correlation function at the crit- 
icality 

It is worth showing how the algebraic decay of the disorder averaged single-particle corre- 
lation function at the critical point arises in d > 1. The Bose operator correlation function 
averaged over disorder is 

~ p^i^e^'t'^i.m-iMm)^ (66) 

where the average is to be taken over the rephcated action. Right at the SF-BG transition, 
the above average becomes 



G{x) ~ po exp -2^K*c* J , (67) 

where rj^ is the anomalous dimension of the phase corelator, defined as 

dlnZK , . 

where the renormalization factor is defined as K{s) = KZk- Since under the scale 
transformation A^A/s, k ^ k/s, it follows that 

dlnZK . s 

^^^-dh^- 

The recursion relation for K may then by definition be written as 

dK 

{2-z-d-rj^)K, (70) 



d\n{s) 



so that at the fixed point we must have r]^ = 2 — z — d. Inserting this into Eq. (67) the 
expression in the exponential becomes 

for X — > oo, and I assumed d — z. The explicit e-dependence in the result arises from the 
integration in > 1. Adhering to the philosophy of minimal substraction we should set 
e = in the last equation, so that the only dependence on dimensionality is the one implicit 
in the value of K* and c* at the fixed point. Thus Eq. (17) obtains. 
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Figure captions: 



Figure 1: The renormalization group flow in d > 1. Close to d = 1 the fixed point which 
controls the transition between the superfluid — > 0, — > 0) and the localized {D —>■ oo, 
K — > oo) phases is at a small value of disorder D* ~ e. Inset a) depicts the flow in d < 1, 
where only the localized phase is stable, and the repulsive flxed point lies in the unphysical 
region D < 0. Precisely a,t d — 1 the flxed point lies on the marginal line ai D — (inset 
b). 
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